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Abstract 

We investigate by molecular dynamics simulations a continuous isotropic core-softened potential 
with attractive well in three dimensions, introduced by Franzese [J. Mol. Liq. 136, 267 (2007)], 
that displays liquid-liquid coexistence with a critical point and water-like density anomaly. Besides 
the thermodynamic anomalies, here we find diffusion and structural anomalies. The anomalies, 
not observed in the discrete version of this model, occur with the same hierarchy that characterizes 
water. We discuss the differences in the anomalous behavior of the continuous and the discrete 
model in the framework of the excess entropy, calculated within the pair correlation approximation. 

PACS numbers: 64.70.Pf, 82.70.Dd, 83.10.Rs, 61.20.Ja 



1 



I. INTRODUCTION 



Experiments for water allow to locate the line of temperatures of maximum density 
(TMD) below which the density decreases with decreasing temperature at constant pressure, 
instead of increasing as in the majority of fluids.^ Experiments for Te,^ Ga, Bi,- S,—'^ and 
GeisTess,^ and simulations for silica,-'^'^^'^ siliconii and BeF2,- show, as well, the same 
density anomaly. 

Water also has dynamic anomalies. Experiments for water show that the diffusion con- 
stant, D, increases on compression at low temperature T up to a maximum D^i^^{T) at 
P = PoiaaxiT). The behavior of normal liquids, with D decreasing on compression, is 
restored in water only at high P, e.g. for P > Pomax ~ 1-1 kbar at 10 °C.- Numerical 
simulations for SPC/E watev^ recover the experimental results and show that the anoma- 
lous behavior of D extends to the metastable liquid phase of water at negative pressure, 
a region that is difficult to access for experiments.— i^^ii^ii^ In this region the diffusivity 
D decreases for decreasing P until it reaches a minimum value D,^ia{T) at some pressure 
-PDmin(^)) a^d the normal behavior, with D increasing for decreasing P, is reestablished 
only for P < Po^i„(T).i^^^^ 

In order to shade some light in the relation between the thermodynamic and dynamic 
anomalies, Errington and Debenedetti investigated the pressure-dependence of the structure 
of SPC/E water by measuring (i) the tendency of pairs of molecules to adopt preferential 
separations, by means of a translational order parameter, and (ii) the extend to which a 
molecule and its four nearest neighbors assume a tetrahedral arrangement, by means of an 
orientational order parameter.— They found a region of the P-T phase diagram, called the 
structural anomalous region, where both order parameters decrease for increasing density, 
i.e. the liquid is more disordered at higher density, in contrast with the behavior of normal 
liquids. They showed that the structural anomalous region encompasses the region where the 
diffusion is anomalous, and that the latter includes the region where the density is anomalous. 
Next, Shell, Debenedetti, and Panagiotopoulos used the two structural order parameters to 
study the case of a computer model for silica. Both order parameters demonstrated to be 
anomalous, but the hierarchy of anomalies for silica is different compared to water. ^ For 
silica the diffusion anomalous region is the wider one, containing the structural anomalous 
region which in turn has inside the density anomaly region. 
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The study of core-softened potentials generates models that are computationally (some- 
times even analytically) tractable and that may retain some qualitative features of net- 
work forming fluids such as water and silica. What they have in common with the water 
potential is the two interaction scales that in the case of water are the van der Waals 
and the hydrogen bond interaction distances. In some cases, theoreticaU^ii^'22i2ii22ii22ii2ii2^ 
and experimental^ results for core-softened potentials show the presence of anomalies. In 
other cases, in particular for the discontinuous shouldered well (DSW) potential, no density 
anomaly is observed.— >2^'2£i22 The DSW potential consists of a hard core, a square repulsive 
shoulder and an attractive square well (dashed line in Fig. [1]). 

In this paper, we investigate the presence of water-like anomalies in a continuous shoul- 
dered well (CSW) potential introduced by Franzese^^ (continuous line in Fig. [T]). Although 
the model is similar to the DSW version, we find density, diffusion, and structural anomalies 
and we observe that they occur with the same hierarchy as in water. In order to get a better 
understanding of the reasons why the CSW potential displays water-like anomalies, while 
the DSW one does not, we apply recent ideas that have emerged on the connection between 
the anomalies and the excess entropy.-'^^'^ 

Since the model studied here resembles effective potentials for complex liquids or liquid 
metals, our results suggest that the water-like hierarchy of anomalies could be experimentally 
found in systems where pressure-induced bond shortening, rather than the change in the 
coordination number, is the origin of the density anomaly.— 

The outline of the paper is as follows. We present the details of the model in Sec. [Tll 
and the details of the simulations in Sec. Illli In Sec. HV] we show and discuss the results 
obtained from simulations and in Sec. |V]we summarize and give the conclusions. 



II. THE MODEL 

The CSW model^ studied here is composed by a system of identical particles interacting 
through the potential (Fig. [1]) defined as 



where Ur and Ua are the energy of the repulsive shoulder and of the attractive well, respec- 
tively, a is the diameter of particles, Rr and Ra are the repulsive average radius and the 
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exp [A (r - Rr) /a] 
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distance of the attractive minimum, respectively, A is a parameter related to the slope of the 
potential at Rr^ 6\ is the variance of the Gaussian centered at Ra- The many parameters of 
this potential allow it to assume several forms, similar to the potential proposed by Cho,— 
ranging from a deep double wells^i^i^ to repulsive ramp-like potentials,— making this 
equation malleable for studying different isotropic models for systems such as colloids or 
liquid metals. Details about the role played by these parameters may be found in Ref.— . 
The set of parameters used here is the same as in Ref.—: Ur/Ua = 2, A = 15, Rn/a = 1.6, 
Ra/cl = 2, ((5yi/a)^ = 0.1, with cutoff = 3.0a. To have a continuous function at the cutoff 
distance, a constant C/Ua = 0.208876 and a linear term Ar/a, with A/a = —0.0673794, 
are added.— The CSW potential may be interpreted as a continuous version of the DSW 
potential investigated in Refs.— The corresponding DSW potential is represented 
by a square shoulder of size w^/a and a well of size wp^/a and depth Ur + Ua illustrated 
by the dashed line in in Fig. [1] (the parameters in figure are wn/a = 0.612, WA/a = 0.767, 
Ur/Ua = 2.2). 

III. DETAILS OF SIMULATIONS 

In order to investigate the interparticle CSW potential Eq. ([1]), we employ a molecu- 
lar dynamics (MD) technique.— We consider 1000 particles in a cubic box with periodic 
boundary condition. We perform simulations in the NVT ensemble with the Nose-Hoover 
thermostat^ and coupling parameter equal to 2. Pressure, temperature, density, and dif- 
fusion are calculated in dimensionless units, P* = Pa^ /Ua, T* = ksT /Ua, p* = pa^, and 
D* = D{m/a'^UAf''^. 

Positions, velocities, and configurational energy of the particles are stored for every 100 
steps, during a total time of 10^ steps, given a total of lO'* stored configurations. Pressure, 
diffusion, and order parameters are calculated over these configurations, discarding the first 
2 X 10^ for equilibration purposes. The time step used is 2 x 10"'^ in units of {c?m/U a)^^'^ (of 
the order of ~ 1.3 x lO^^^s for argon-like atoms and ^ 1.1 x lO^^^s for water-like molecules). 

The model displays^^ a phase diagram with two first-order phase transitions ending in 
critical points: one transition between gas and low-density liquid (LDL) phases and another 
between LDL and high-density liquid (HDL) phase. The LDL-HDL phase transition is 
metastable with respect to the crystal phase. 
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For temperatures and densities inside the metastable region the positions, velocities, and 
configuration energy are stored for every 10 steps and the quantities are calculated with 
the configurations stored until the drop in configuration energy and virial, discarding the 
equilibration time. The time before the drop in energy and virial characterizes the lifetime 
of the metastable phase. 

IV. RESULTS AND DISCUSSION 

A. Density anomaly 

In the vicinity of the LDL-HDL critical point the system shows density anomaly.— The 
isochores have a temperature of minimum pressure in the interval 0.185 < p* < 0.215, 
corresponding to the temperature of maximum density (TMD) at constant P [Fig. [2]^a)]. 
Inside the TMD line the density decreases as the system is cooled at constant P. For densities 
in the range 0.2 < p* < 0.215, the TMD line is negatively sloped, while it has positive slope 
for 0.185 < p* < 0.2. Our results for the TMD line in the ranges 0.55 < T* < 0.61 and 
0.25 < P* < 0.30 [Fig. El^d)] compares well, qualitatively, with the experimental data for 
water,- showing a negative slope [Fig. Et^b)]. The data in all the supercritical region [Fig. 
I21^a)] resemble the simulation results for SPC/E water, with a retracing TMD line [Fig. 
|2]^c)]. Nevertheless, in the present case the TMD line is at P > for all its length, while for 
SPC/E water it extends to negative pressures. 

B. Diffusion anomaly 

The diffusion coefficient is calculated using the mean-square displacement averaged over 
different initial times. 



(Ar(t)^) = ([r(to + t)-r(to)]^). 



(2) 



From Eq. ([2]), the diffusion coefficient may be obtained as follows: 



t— >oo 



(3) 
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While for normal fluids diffusivity decreases monotonically with increasing density at 
constant temperature, for the model Eq. ([1]) this is the case only for temperatures T* > 0.64 
(Fig. [3]). For T* < 0.64, the diffusion coefficient has three regions: (i) For p < pDmin, 
D decreases as p increases as expected for normal fluids, (ii) For p Drain < P < PDmax, 
D increases with density - an anomalous behavior. This region is called the anomalous 
diffusion region, (iii) For p > pomax, the normal behavior is restored with D decreasing as 
density increases. 

The line of PDmni^iT) in the P-T phase diagram [Fig. EJ^d)] is consistent with the diffusiv- 
ity maxima observed in experiments for water [Fig. EJ^b)]. Moreover, both lines of diffusion 
extrema, PDmaxiT) and PDminiT) [Fig. [2](a)], resemble the results for detailed models of 
water— li^ii^ [Fig. [2|^c)], silica,-ii^i^ and other isotropic potentials.— We flnd that our 
Dmin, like other isotropic models,— i^^i^^i^^ occur at positive P. 

C. Structural anomaly 

The two quantities used for studying the structural behavior of the system of particles 
interacting through the potential in Eq. ([1]) are the translational order parameter, t, and 
the orientational order parameter, Qq. The translational order parameter is deflned as^>i^i^ 



where ^ = rp^^^ is the distance r in units of the mean interparticle separation p^^^'^, is the 
cutoff distance set to half of the simulation box times p~^^^, as in Ref."^-, g{^) is the radial 
distribution function proportional to the probability of flnding a particle at a (reduced) 
distance ^ from a reference particle. For an ideal gas g = 1 and t = 0. In the crystal phase 
g ^ 1 over long distances and t is large. 

For normal fluids, t increases with increasing density. We flnd this monotonic behavior 
only for T* > 1.8 (Fig. Hj). For T* < 1.8 we observe that t has a maximum at ptmax{T), 
decreases for increasing p > ptmaxiT), reaches a minimum at ptiamiT), and recovers the 
normal increasing behavior for p > ptmin{T) (Fig. |4]). 




(4) 



For each particle i = 1, . . . , we calculate^'^^'^'^'^'^'^ 
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with i = 6. The quantity {Y^^^, = 




harmonics Y^m calculated over the vectors rij{9ij, with j = 1, . . . , k, connecting particle 
i with its k nearest neighbors j. We use = 12 as in Ref.^^. The local orientational order— 
of the whole system is calculated as the average over the particles, 



Our study reveals an anomalous behavior also for the orientational order parameter, 
Qq. For normal fluid, Qq increases monotonically with increasing density. Instead, for the 
present model we find that Qq is non monotonic, displaying a maximum at pgmax; for all 
the considered temperatures (Fig. [5]). 

Since pgmax lies between the densities which bound the extrema in t, ptmax and ptmin, we 
call structural anomalous region the interval between pg^aax and ptmm- In this region both 
parameters t and Qq are anomalous, and the liquid becomes less ordered with increasing 
density, in contrast with the behavior of normal liquids. 

D. The hierarchy of anomaUes and order map 

The relation between the several anomalies presented for the CSW potential is illustrated 
in Fig. El The TMD line lies between the diffusivity extrema (DE) lines, that are included 
within the structural anomalous region bounded by the curves of maxima of Qq and minima 
of t. The hierarchy of anomalies found here is the same reported for the SPC/E water— and 
other two-scales potentials.— 

The order map in the t — Q^ plane (Fig. [7]) resembles the one observed for SPC/E water,— 
silica,- and other two-scale potentials2ii2i since it has an inaccessible region. Differently from 
water, and similar to other two-scale potentials,—"^ the translational and orientational order 
parameters are not coupled into the structural anomaly region. Indeed, for the densities 
belonging to the structural anomalous region, t and map onto a two dimensional region 
(Fig. [7]), in opposition to the water case, in which they map onto a single line.— 




(6) 



i=l 



For a crystal, Qq is large, while Qjf = l/\/k for an ideal gas. 
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E. The Widom line and the minima in t 



The so called Widom line is defined as the locus of the maximum of response functions 
close to the critical point at the fluid phase. Here we calculate the isothermal compressibility 



As expected, the Widom line (diamond symbols in Fig. [8]) is a continuation of the liquid- 
liquid coexistence line (crosses in Fig. [8]). Our data show that the Widom line coincides with 
the minima in the translational order parameter t, at least in the vicinity of the liquid-liquid 
critical point. This is consistent with the observation that compressibility, structure factor, 
g{r) and the parameter t are related. In particular, we find that kt has its maximum where 
the translational order is minimum. Hence, the largest variation of volume at constant T 
for increasing P occurs at the pressure where the behavior of the structural parameter t 
changes from anomalous to normal. To our knowledge these are the first data showing this 
relation between kt and t. 

F. Excess entropy and anomalies 

Why the DSW potential has no water-like anomalies and its continuous counterpart, the 
CSW potential, does? We can gain some understanding by analyzing the density dependence 
of the excess entropy. It has been shown® that the behavior of the excess entropy as a function 
of density is linked to anomalies in density, diffusion and structure. 

Extrema in {dS/dp) correspond to an thermal expansion coefficient equal to zero and 
therefore to density extrema.- The thermodynamic condition that gives rise to the density 
anomaly can be written as Sex = {ds'^^/dln p)j, > 1, as shown by Errington et al. in Ref.— . 
Here s'^^ is the excess entropy, defined as s^^ = s — s^^, i.e. the difference between the 
entropy s of a real fluid and the entropy s^^ of an ideal gas at the same T and p, due to the 
correlations between the position of the particles of the real fluid. 

Based on the proposition of Rosenfeld^ that the logarithm of the diffusion coefficient 
D is proportional to the excess entropy, a non-monotonic behavior in the excess en- 
tropy would imply a non-monotonic behavior in D. Following the empirical Rosenfeld's 
parameterization,-^ Errington et al^'^ have also observed that the diffusion anomaly can be 



(Fig. ED, 




(7) 
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predicted by Sex > 0.42. 

Finally, they argue that Sex > is a good estimate for determining the region where struc- 
tural anomaly occurs, because for normal fluids the excess entropy decreases for increasing 
density at constant temperature.— Recent works have explored the new possibilities one 
can achieve through the connections between excess entropy and structure, not only in the 
matter of isotropic fiuids^'^'^ but also in water,— silica,- and BeF2-—^ 

To calculate the excess entropy, one should count all the accessible configurations for a 
real fluid and compare with the ideal gas entropy. This calculation is not straightforward 
and can approximated by 



since S2 is the dominant contribution to excess entropy^i^i^ and it is proved to be between 
85% and 95% of the total excess entropy in Lennard- Jones systems.—*^ The two-body 
contribution S2 depends only on the radial distribution function g{r) and the density, giving 
a direct connection between structure and thermodynamics. The excess entropy and the 
translational order parameter are linked because both depend on the deviation of g{r) from 
unity. The relation between excess entropy and the orientational order parameter, however, 
depend on the symmetries of the structures in the high and low density limits. Indeed, it 
turns out that the extent to which the pair correlation entropy is sensitive to orientational 
order determines the behavior of the structurally anomalous region.- 

However, as shown in Ref.— S2 overestimates the excess entropy, s'^^ < S2, and Sex < 
S2 = {ds2/d\n p)rp. Therefore, the inequalities S2 > 0, S2 > 0.42, and S2 > 1, overestimate 
the region of structural, diffusion and density anomalies, respectively. Nevertheless, S2 has 
been shown to give estimates of the anomalous regions in qualitative agreement with the 
estimates based on Sgx.— 

In order to compute S2 for the DSW potential, we simulate 500 particles inside a cubic 
box with periodic boundary conditions. The particles interact through the DSW potential 
(dashed line in Fig. [1]), with parameters wr/o = 0.612, wa/o = 0.767, and Ur/Ua = 2.2. 
The equilibration and production times are 350 and 650, respectively, in reduced units. 
To achieve the desired temperature, we rescale the velocities until the equilibration. At 
equilibrium, we simulate the system in the NVE ensemble. 

From the analysis of S2 and S2 for the DSW potential [Fig. [9](a) and Fig. [9]^b)], we 
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observe that S2 has a maximum, not monotonic with T, that is always smaller than 1. 
Therefore, Sex < 1 and no density anomaly is expected, in agreement with Franzese et al..— 

We find densities where the DSW potential has S2 > 0.42, for 0.60 < T* < 0.75, and 
S2 > 0, for all the four temperatures studied here [Fig. El^b)]. These values apparently 
suggest the presence of diffusion and structural anomaly, respectively. However, since S2 > 
Sex and the Rosenfeld criteria holds with 30% uncertainty,'^ our results for S2 in Fig. Mjo) 
give no final answer about the presence of diffusion anomaly. A more detailed study, beyond 
the goal of this work, on diffusion and structural anomaly for the DSW potential would be 
necessary to clarify these points. 

Our results on S2 and S2 for the CSW potential are consistent with our analysis of the 
anomalies described in the previous sections [Fig. M^c) and Fig. El^d)]. We observe that 
S2 displays a maximum that increases for decreasing temperature. For the four studied 
temperatures, we find that S2 > 0.42 (hence, greater than zero), suggesting the occurrence 
of both diffusion and structural anomalies. Since we find a range of densities where S2 > 1 
only for T* < 0.75, the result suggests the presence of density anomaly only below this 
temperature. This is qualitatively consistent with the water-like hierarchy of anomalies that 
we have reported in the previous sections. However, the comparison between Fig. [9](d) and 
Fig. [6] shows that the criteria based on S2 overestimate the temperatures and densities 
where the anomalies appear. We interpret this discrepancy as the effect of using S2, instead 
of Sex, in the criteria for the anomalies, with S2 > Sex- 

Therefore, the excess entropy analysis is a useful tool to distinguish between the DSW and 
the CSW potential. In the discontinuous case the increase of the excess entropy with density 
at constant T is not enough to give rise to the density anomaly. In the continuous case, 
instead, its increase is large enough at low T. For the diffusion and structural anomalies, 
the analysis based on the approximate expression S2 is less clear. However, the comparison 
with the direct calculations of the anomalies in the CSW case suggests that only for the 
CSW potential the increase of excess entropy is enough to reach the regimes with diffusion 
and structural anomalies at the considered temperatures. 

From Fig. [9] one observes that for the CSW potential S2 is more sensitive to density 
variations than for the DSW potential when p* < 0.24, i. e. when the average interparticle 
distance r/a= {p*)~^^^ > 1.612. This is consistent with the fact that the DSW potential is 
constant for 1.612 < r/a < 2.379, corresponding to average density p* < 0.24 and p* > 0.07, 
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respectively. Therefore, for the DSW potential in this regime of density one can expect 
a small variation of g{r) and, as a consequence, of S2. Instead, the CSW potential for 
r/a > 1.612 (p* < 0.24) has a large variation, equal to 100% of its attractive energy depth, 
implying a sensible variation in g{r) and S2- 

It is worth noticing that, recently, Netz et alM have argued that, in potentials with two 
characteristic scales, the region of density anomaly can be observed only if it appears at a 
temperature such that ksT is larger than the discontinuity in the interaction energy (e.g., 
> Ua + Ur in the DSW in Fig. [1]). This is consistent with our results for the CSW 
potential, showing the density anomaly at ksT = 0.65 Ua > 0, because the CSW potential 
has no discontinuity in the interaction energy. On the other hand, for the DSW potential 
with discontinuity Ua + Ur = 3.2 Ua, since no density anomaly is found for ksT > 3.2 Ua, 
this criterion excludes that the anomaly could be found at lower temperatures. 

V. SUMMARY AND CONCLUSIONS 

We studied a three dimensional system of particles interacting through a continuous 
isotropic interparticle pair potential recently proposed by Franzese.— The potential has a 
hard-core, a repulsive shoulder and an attractive well. This continuous shouldered well 
(CSW) potential resembles the discontinuous shouldered well (DSW) potential studied in 
Ref.— and shows new features that are absent in the DSW potential. 



TABLE I: Critical temperatures and T^^, pressures P^^ and Pq^, and densities Pq_^ and p^^, 
for the gas- liquid critical point Ci and the liquid-liquid critical point C2, for the CSW potential. 
The quantities are expressed in dimensionless units. 



Critical point 




p* 


P* 


Gas-Liquid Ci 


0.95 ± 0.06 


0.019 ±0.008 


0.08 ± 0.03 


Liquid-Liquid C2 


0.49 ±0.01 


0.285 ± 0.007 


0.247 ± 0.008 



As its discontinuous counterpart, the CSW potential displays a phase diagram with a gas- 
liquid phase transition ending in a critical point, and a LDL-HDL phase transition ending 
in a liquid-liquid critical point. Table [T] shows the values of the critical parameters for the 
two critical points. The LDL-HDL phase transition is metastable with respect to the crystal 
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phase.— In the P-T phase diagram the LDL-HDL phase transition hne has positive slope, as 
in the discontinuous model,— and in the ramp-potential.— This feature suggests that these 
models are describing systems different from water, because for water and water-like models 
the slope of the liquid-liquid phase transition is expected to be negative.— i^i^^^ 

In contrast with its discontinuous counterpart, the CSW potential has density anomaly 
(FiglH]).— The TMD line bends toward the LDL-HDL critical point at high P, as also seen 
for other isotropic potentials.— 

We find that the CSW potential has also diffusion and structural anomalies (Fig. [8]), 
as water and silica. Since water— and silicar' have different order of these anomalies, we 
investigate the cascade of anomalies for the CSW potential. We find that for the CSW 
model the structural anomalous region encompasses the diffusion anomaly region, inside 
which the density anomaly region is observed, as in water. This is consistent with previous 
analysis of other two-scale isotropic potentials^ii^^. 

The loci of structural anomalies and the two phase transition lines in the P-T phase 
diagram display a positive slope. The TMD line has positive slope at low P and negative 
slope at high P. By decreasing T, we find that three quantities - the minima of the structural 
order parameter, the Widom line, and the maxima of diffusivity - approach the LDL-HDL 
critical point value, suggesting that the high-T behavior of these quantities could give an 
indirect indication of the location of the LDL-HDL critical point in real systems, if it is 
present. In particular, close to the LDL-HDL critical point, the Widom line coincides with 
the minima of the translational structural order parameter. 

To understand the difference between the DSW and the CSW potential we perform an 
excess entropy analysis, approximating the excess entropy by the pair correlation entropy 
S2- Comparison of this analysis with the direct calculation of the anomalies suggests that 
the approximate excess entropy S2 satisfies criteria with empirical thresholds higher then 
those predicted for the exact excess entropy s^^. 

Interestingly, the excess entropy analysis for the DSW and the CSW potential emphasizes 
a relevant difference for the appearance of the anomalies: the DSW potential is constant 
in the density regime of the anomalies (approximately for 0.17 < p* < 0.25, i. e. average 
interparticle distance 1.6 < r/a < 1.8), while the CSW potential has r-dependent soft re- 
pulsion. This observation could clarify also why potentials, such as the ramp in Ref.— , show 
water-like anomalies, because they have appreciable distance dependence (soft repulsion) for 
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interparticle separation between the repulsive and attractive length. 

Our work confirms that water-like anomalies can be present also in systems that, unlike 
water, have no directional bonds, as discussed in previous works.— >i2ii22i2ii2Ii2^iMi2^ The CSW 
potential, and other two-scales isotropic models^'i^'^i^'^i^i^'^i^, have features that are 
different from water, such as the positive slope in the P-T phase diagram of the LDL- 
HDL phase transition, or the locus of diffusion minima at positive pressure, instead of 
negative.— These differences emphasize the possible existence of liquids with anomalies 
induced by bond shortening at high P, rather than the change in the coordination number 
as in water or silica. This conclusion is consistent with recent results on metallic glasses 
showing polyamorphism.— 

We conclude by remarking that, although a similar cascade of anomalies was observed in 
other two-scales potentials,— '2^'^ the CSW potential is the first with a strong shoulder that 
displays these anomalies. This feature make the CSW potential suitable to study liquid 
metals, such as Cs, Mg, Sr, Ba,— and colloids,— where effective interactions with large 
shoulder are derived in some cases. By varying the parameters of Eq. ([T]), the potential 
changes from deep double wells to a repulsive ramp, and by varying the power of the core of 
Eq- dl]) ^ here equal to 24 to mimic a hard-core - the stiffness of the core changes, making 
the Eq. ([1]) an interesting functional form to model and interpolate pressure-dependent 
effective interactions. 
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Captions to the figures 



Fig. 1: (Color online) Interaction potentials studied in this work. Dashed line represents 
the discontinuous shouldered well (DSW) potential.— '^^I'^Si'^ Continuous line represents the 
continuous shouldered well (CSW) potential introduced in Ref.— . The parameters are ex- 
plained in the text. 

Fig. 2: (Color online) (a) Pressure-temperature (-P-T) diagram for the CSW model. 
Lines correspond to isochores with, from bottom to top, p* = 0.175, 0.18, 0.185, 0.19, 0.2, 
0.21, 0.215, 0.22, 0.23, 0.24, 0.25, 0.26, 0.27, 0.28, 0.29, 0.3, 0.31, and 0.32. Triangles 
represent spinodal lines for the LDL-HDL phase transition, with LDL at low P and HDL 
at high P. We estimate the LDL-HDL critical point where the spinodal lines converge 
(large filled circle). Line with crosses represents our estimate of the liquid-liquid coexistence 
line. The TMD (bold continuous) line bends toward the LDL-HDL critical point at high P. 
Dashed lines bound the region where the diffusion anomaly occurs (see section |IVB|) . (b) 
Experimental data for water anomalies adapted from Angell et al.r- Circles denote the line of 
temperatures of maximum density (TMD) at constant P. Squares mark where the diffusion 
has a maximum value with increasing P at constant T, PDmax(^)- (c) Simulation data 
for SPC/E water adapted from Netz et al..— Squares mark Pomax where the diffusion has a 
maximum value, -Dmax? at constant T, and diamonds mark Pomm ioT the local minima, -Dmin- 
Circles locate the TMD line, (d) Zoomed region from panel (a), showing good qualitative 
agreement between our simulations and the experiments. 

Fig. 3: (Color online) The diffusion coefficient against the density for several isotherms. 
For the range of densities bracketed within the dashed lines, the particles move faster under 
compression for temperatures lower than 0.64. This is the opposite behavior which one 
expects for normal fluids. Dashed lines are guides for the eyes connecting p_D„.„ and pomax- 

Fig. 4: (Color online) The translational order parameter as a function of density. While 
for normal fluids compression leads to increase the translational order parameter, for the 
model of Eq. ([T]) this is the case only for high temperatures (T* > 1.8). Dashed lines bound 
the region where t behaves anomalously. 

Fig. 5: (Color online) The orientational order parameter against density. We observe 
that Qq has a maximum at pQmax, meaning that decreases under compression for some 
range of densities. The Qq maxima lie between the extrema points of the translational order 
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parameter, t. Dashed line marks the location of maximum Qq. 

Fig. 6: (Color online) Temperature- density plane containing all the anomalies found for 
the CSW potential. The TMD line bounds the innermost region with the density anomaly 
behavior. This region is surrounded by the D extrema lines, that encompass the region 
with diffusion anomaly. The out-most anomalous region, including the first two, is between 
curves B and A, where the system exhibits an anomalous behavior in structure as shown 
by the order parameters t and Qq. The curve C marks the maxima in t occurring where Qe 
has a normal behavior. 

Fig. 7: (Color online) The t — Qq plane or order map. The arrows indicate the direction 
of increasing density. Each line correspond to an isotherm and from top to bottom they 
are T* = 0.49, 0.55, 0.60, 0.62, 0.65, 0.70, 0.75, 0.80, 0.90, 1.0, 1.3, 1.5, 1.7, and 1.8. By 
increasing the density, at low p both order parameters increase (normal fluid behavior), 
then at intermediate p they both decrease (structural anomaly region), then at higher p 
the orientational Qq decreases, while the translational t increases. As in the case of SPC/E 
water, silica and other two-scales potentials the region with high Qq and low t is inaccessible. 
The inaccessible region is limited by a straight line tmin = a + bQQ with a = —2.86 ± 0.02 
and b = 12.8 ± 0.1. 

Fig. 8: (Color online) Pressure-temperature phase diagram merging all the results found 
for the CSW potential. The meaning of the lines is described in the legend, where DE stands 
for diffusivity extrema, LL for liquid-liquid, and LG for liquid-gas. See the text for more 
details. 

Fig. 9: (Color online) (a) Pair contribution of excess entropy, S2, for the DSW potential 
(dashed line in Fig. [1]) against density at constant T. Circles are simulated data and lines 
are fifth order polynomial fit from data, (b) S2 = {832/ din p)rp is shown for DSW potential. 
Panels (c) and (d) show the results for the CSW model. Horizontal lines mark the threshold 
value for anomaly in density p, diffusion Z), and structure, as explained in the text. Solid, 
dotted, dashed, and dotted-dashed hues correspond to temperatures T* = 0.55, 0.60, 0.75, 
and 0.90 in all panels. 
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